ABSTRACT 



We study the electronic states of isolated fullerene anions Cqq (1 < n < 6) 
taking into account the effective interaction between electrons due to exchange of 
intramolecular phono ns. If the vibronic coupling is strong enough such an effect may 
overwhelm Hund's rule and lead to an ordering of levels that can be interpreted as 
on-ball pairing, in a manner similar to the pairing in atomic nuclei. We suggest 
that such effects may be sought in solutions of fulleride ions and discuss recent 
experimental results. 
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I. INTRODUCTION 

The discovery of superconductivity in alkali-metal-doped fullerenes 1 K 3 C 60 and 
RD3C60 has raised interesting questions about the electron-phonon coupling in such com- 
pounds and its interplay with Coulomb repulsion. Ceo is a highly symmetrical molecule 
i.e. it is a truncated icosahedron and its electronic lowest unoccupied molecular orbitals 
(LUMO) are threefold degenerate 2 ' 3 ' 4 . They form a T lu representation of the icosahedral 
group Ih- Filling the LUMO in Cqq anions leads in a naive picture to narrow, partially 
filled bands in the bulk fullerides. The bandwidth W is determined by the hopping be- 
tween the C 60 molecules which are quite far apart and W 0.5 eV. The coupling of some 
K g phonons with electrons residing in the Ti u orbital has been suggested to be responsi- 
ble for the superconductivity 5 ' 6 ' 7 . The Coulomb repulsion also may be important on the 
ball 8 . Several authors 9-12 have undertaken the study of the Jahn- Teller distortion that is 
expected in the fullerene anions. In such calculations one considers the electrons as fast 
degrees of freedom and the phonon normal coordinates are treated as static 13 . 

In this paper we investigate the interplay between the electronic and phononic degrees 
of freedom on an isolated fullerene anion. We study an effect that goes beyond the Born- 
Oppenheimer approximation which is the modification of electronic levels due to phonon 
exchange. We obtain the lifting of degeneracy by a perturbation calculation in the case 
of an undistorted anion. The ordering of levels can be described as "anti-Hund"' rule. 
Our calculation is very close in spirit to the standard treatment of the electron-phonon 
coupling in superconducting metals. Here we argue that the energy scale of this effect 
may be comparable to that of the Jahn- Teller effect. This is because the phonons have 
high frequencies as well as medium to strong coupling to electrons. The effect we observe 
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may be sought by spectroscopy of solutions of fullerides in liquid ammonia, for example. 
We discuss the opposite effect of Coulomb interaction, leading to Hund's rule in ordinary 
situations. Finally we point out that experimentally observed spectra may be at least 
partially explained by our crude calculation. 

In section II we discuss the electron-phonon coupling, section III gives our results 
for the electronic levels of the anions, section IV discuss the competition with Coulomb 
repulsion and section V contains a brief discussion of near-IR spectroscopic measurements 
and EPR experiments on solutions of fulleride anions as well as our conclusions. 

II. THE ON-BALL ELECTRON-PHONON INTERACTION 

The electronic structure of % electrons in the Ceo molecule is well known to be given 
by a simple Hiickel calculation. The levels are labeled 4 by the irreducible representations 
(irreps) of the icosahedron group Ih- One important property has to be noted: three of 
the Ih irreps are the /=0,1,2 spherical harmonics of SO(3) which do not split under the 
Ih group. They are commonly named A 5 , Ti„, H g . In addition there is also the twofold 
spin degeneracy. 

In the ground-state of the neutral Ceo molecule all levels up to H u included are 
completely filled thus building a singlet state |*o)- The LUMO are the six T lu states. 
These are occupied upon doping with extra electrons and the ground-state becomes then 
degenerate. One then expects the Jahn- Teller effect to distort the anion and lift this orbital 
degeneracy 13 ' 14 . We focus on another effect which goes beyond the Born-Oppenheimer 
approximation in the sense that nuclear motions are crucial for its very existence: the 
coupling of the T lu electrons to the vibrational modes of the molecule (also referred to 
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as phonons). Phonon exchange between electrons leads to an effective electron-electron 
interaction that competes with Coulomb repulsion and may lead to anti-Hund ordering of 
energy-levels. 

For simplicity we treat this effect assuming the absence of the Jahn- Teller distortion. 
The next logical step would be to compute first the distortion pattern of the anion under 
consideration, then obtain its vibrational spectrum and the electron-phonon coupling in 
the distorted structure and then compute again the effective electron-electron interaction. 
As a first investigation of electron-phonon coupling we use a perturbation scheme suited to 
degenerate levels we will derive an effective electron-electron interaction with the assump- 
tion that filled states lying below the T lu level remain frozen so that intermediate states 
involve only T lu -T lu excitations. Indeed the H u -T lu gap is 2eV whereas maximum 
phonon energies are ~ 0.2eV. 

A typical electron-phonon interaction term reads: 



Here X a are normal coordinates, the subscript referring both to the irrep and to the row 



the Tiu (^=1) level, mi taking one of the m=— 1,0,1 values, and / Qmim2 are complex 
coefficients. The cj^ operators transform as 1=1 \l,m) vectors under Ih symmetries, 
and their conjugates c mcr transform as (— l) m+1 |Z, — m) vectors. The (— l) m2+1 cJ nicr c_ m2Cr 
products transform then as members of the Ti u x Ti u representation, which in the Ih 
group splits as: 




(H7!lTO2 « "ll" "12^' 



a,mi,rri2,cr 



in the irrep they belong to, cj, 



is the creation operator for an electron with spin a in 



Tlu X Tlu = A ff + Tig + Hg. 
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This selects the possible vibrational modes T lu electrons can couple to. In fact, only H g 
modes split the degeneracy 7 . 

Let us consider a particular fivefold degenerate multiplet of H s modes. Their normal 
coordinates will be labelled X m , m ranging from -2 to +2. Since H g appears only once 
in the product Ti u x Ti u , the interaction is determined up to one coupling constant g by 
the usual formula for the coupling of two equal angular momenta to zero total angular 
momentum: 

W = gJ2(- 1 ) mX m$-m (1). 

m 

The X m may be chosen such that X^ = (— l) m X_ m and have the following expression in 
terms of phonon operators: 

X m = -L (a m + (-l) m al m ) (2) 

whereas the $ m are the irreducible 1=2 tensor operators built from the c^c products ac- 
cording to: 

d> m = £(1, 1, 2\ mi , m - m lt m){-lf"^cl ia c_ m+mia (3) 

mi 

where (Zi, I2, l\mi, mi-, m) are Clebsch-Gordan coefficients. 

We now consider a doped CgQ~, molecule, < n < 6. Its unperturbed degenerate 
ground-states consist of |^o) to which n T\ u electrons have been added times a zero- 
phonon state. They span a subspace denoted by So. In £ the unperturbed Hamiltonian 
Hq reads: 

m,a m 
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where e tlu is the energy of the T lu level, hu is the phonon energy of the H 3 multiplet 
under consideration. Within So the effective Hamiltonian up to second order perturbation 
theory is given by: 

H eff = E Q P + PoWPo + P W(1 - Pq ) 1 - (1 - P )WP 

-t^o — Hq 

where Po is the projector onto So, E is the unperturbed energy in this subspace which is 
just the number of doping electrons times et lu - The linear term in W gives no contribution. 
Using expressions (1) and (2) for W and X m one finds: 

H eff = H °-^i E (-^ma^-ma 2 (4) 
m,cr 1 <J2 

where we have now included spin indices. We can now use equation (3) to express H e ff 
as a function of c and operators and put it in normal ordered form using fermion 
anticommutation rules. In this process there appears a one-body interaction term which 
is a self-energy term. We will henceforth omit the H term which is a constant at fixed 
number of doping electrons. 

Let us now define pair creation operators A^ which when operating on the vacuum 
|0) create pair states of T iu electrons that are eigenfunctions of L,S, L Z ,S Z , where L,S are 
total angular momentum and spin, and L Z ,S Z their ^-projections. / and s can take the 
values 0,1,2 and 0,1 respectively. This holds also if |0) is taken to be the singlet state |^o)- 

A im = E ( 1 ' 1 ' / l mi ' m_mi ' m )(^' s l (7l ' (T ~ (Tl ' (T ) c ^i^i c ^-^i^-^ ( 5 ) 

mi,<7i 

The quantity A\^J is non-zero only if (/ + s) is even and the norm of A^ \0) is then equal 
to y/2. The inverse formula expressing products as A^ operators is: 

c L 1 a 1 c L 2 a 2 = X^ 1 ' X ' l \ mi > m2 ' mi + m2 )(^' \ a2 ' ai + a ^ A Z\ + +2j ■ ( 6 ) 



l,s 
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As H e f / is a scalar, its two-body part may be written as a linear combination of diagonal 
A\ m ^ A\ m products whose coefficients depend only on / and s: 

J2f(i, s )AZ'az 

ls,ma 

The F{1, s) coefficients are calculated using expressions (4), (3), (6). We then get H e ff in 
final form: 

^ .oot.oo 1 V/l 1 ' Tt 4 k x 1 V4 oot 4 00 I (7\ 

IX -+- ^1 00 ^1 00 - - s±lm ^lm + YO 2-^ 2m 2m I ' ^ ' 

m,a m J 

In this formula N is the electron number operator for the Ti u level; the N term ap- 
pears when bringing H e ff of expression (4) in normal ordered form. In our Hamiltonian 
formulation the effective interaction is instantaneous. 

There are actually eight H s multiplets in the vibrational spectrum of the Cgo molecule. 
To take all of them into account we only have to add up their respective coefficients 
5g 2 /6hu : their sum will be called A. 

III-THE ELECTRONIC STATES OF FULLERENE ANIONS 

We shall now, for each value of n between 1 and 6, find the n-particle states and 
diagonalize H e ff. The Hamiltonian to be diagonalized is that of equation (7) where the 
prefactor is replaced by —A. The invariance group of H e ff is Ih x SU(2). The n-particle 
states may be chosen to be eigenstates of L, S, L 2 , S z and we shall label the multiplets by 
(Z,s) couples, in standard spectroscopic notation ( 2s+1 L stands for (l,s)). The pair (/, s) 
label SO(3) x SU(2) irreps which, as previously mentioned, remain irreducible under 
Ih x SU(2) as long as / doesn't exceed 2; for larger values of / SO(3) irreps split under 
Ih- Fortunately enough, the relevant values of I never exceed 2. Moreover given any value 

7 



Heff = ~6t^{ 



of n, (/, s) multiplets appear at most once so that the energies are straightforwardly found 
by taking the expectation value of the Hamiltonian in one of the multiplet states. The 
degeneracies of the levels will then be (2/ + 1) (2s + 1) . We now proceed to the construction 
of the states. 

• n=l: There are six degenerate 2 P states c^J^o) whose energy is —A. 

• n=2: There are 15 states, generated by applying A\^J operators on \^q). There is 
one 1 S state, nine 3 P states and five 1 D states. 

• n=3: There are 20 states. States of given l,m,s,a can be built by taking linear 
combinations of A^c^^/o) states according to: 

(h, 1, Z|mi, m - m 1 ,m)(s 1 , -, s|tri, a - cr u a)A s ^ c ] m _ miCJ _ ai |* >- 

These states belong to the following multiplets: 4 S, 2 P, 2 D. 

• n=4: There are 15 states, which are obtained by applying A\^ operators on 

/loot 

Ax) l w o/- 

• n=5: There are six 2 P states which are c^^Aqq^Aqq^ |^ ) an d whose energy is —A. 

• n=6: There is one 1 S state whose energy is 0. 

The corresponding energies are given explicitly in Table I and displayed in Fig. 1. 
It is interesting to note that the above treatment of electron-phonon interaction parallels 
that of pairing forces in atomic nuclei 15 ' 16 . Of course in the case of finite fermionic systems 
there is no breakdown of electron number but there are well-known "odd-even" effects that 
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appear in the spectrum. In our case pairing shows up in the 1 S ground state for Cq rather 



than 3 P as would be preferred by Coulomb repulsion i.e. Hund's rule. The construction of 
the states above is that of the seniority scheme in nuclear physics 16 . We note that similar 
ideas have been put forward by V. Kresin some time ago, also in a molecular context 17 . 
The effective interaction that he considered was induced by a core polarization. 

IV-THE EFFECT OF COULOMB REPULSION 

We now consider the Coulomb electron-electron interaction and assume it to be small 
enough so that it may be treated in perturbation theory. To get some feeling of the order of 
magnitude of this repulsion we use the limiting case of on-site interaction i.e. the Hubbard 
model. This Hamiltonian is not specially realistic but should contain some of the Hund's 
rule physics. The two-body interaction now reads: 



where the i subscript now labels the ir orbitals on the Cqq molecule. The quantity U is 
~ 2-3 eV from quantum chemistry calculations 18 Since level degeneracies are split at first 
order in perturbation theory we confine our calculation to this order and have thus to 
diagonalize the perturbation within the same subspace Sq as before. In this subspace it 
reads: 



where greek indices label one-particle states belonging either to |*o) or to the T lu level. 
Let us review the different parts of Wh- Note that since the \^o) singlet remains frozen 
we have the identity: c^Cg = 5 a p if a, f3 label states belonging to l^o)- 




W H = U (a\i)(P\i)m(i\S) ci T 4c 7i c 5T 
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A part involving states belonging to |^ ) only: 



i,a(3 



a, (3 belong to |^o)- This term is thus diagonal within £ and merely shifts the total 
energy by a constant that does not depend on the number of doping electrons. It won't 
be considered in the following. 

-A part involving both states belonging to |*o) and to the Ti u level: 



The sum over (3 is just the density on site i for a given spin direction of all states belonging 
to |^o) which is built out of completely filled irreps. As a result this density is uniform 
and since |^ ) contains 30 electrons for each spin direction it is equal to 1/2. Wh 2 then 
becomes diagonal and reads: 



a, a 

Its contribution is thus proportional to the number of T lu electrons. It represents the 
interaction of the latter with those of the singlet and we won't consider it in the following. 

-A part involving only states belonging to the T lu level: 
Wh 3 has the same form as Wr with all indices now belonging to the Ti u level. Whereas the 
interaction has a simple expression in the basis of \i) states, we need its matrix elements in 



W H2 =U (a\i){i\6)\{Pi\)\ 2 cLcs* cl-rCp-a 



where a, 5 belong to the 



Tiu level whereas (3 belongs to |^o)- It reduces to: 
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the basis of the T lu states. There are in fact two T lu triplets in the one-particle spectrum 
of the Ceo molecule, the one under consideration having higher energy. To construct the 
latter we have first constructed two independent sets of states which transform as x,y,z 
under Ih- These are given by: 

\a) = ^^e^.fl \i) and \a)' = ^^e^.ki \i) 

i i 

where e a are three orthonormal vectors, i labels sites on the molecule, the r, are the vectors 
joining the center of the molecule to the sites while the ki join the centre of the pentagonal 
face of the molecule the site i belongs to to the site i. We assume that the bonds all have 
the same length. These states span the space of the two Ti u triplets. The diagonalization 
of the tight-binding Hamiltonian in the subspace of these six vectors yields then the right 
linear combination of the \a) and \a)' states for the upper lying triplet. From the x,y,z 
states one constructs 1=1 spherical harmonics. We then get the matrix elements of Wh 3 
in the basis of Ti u states. As So is invariant under Ih operations and spin rotations, Wh 3 
which is the restriction of Wh to So is invariant too. It may thus be expressed using the 
A, A^ operators by using formula (5) in the same way as the phonon-driven interaction 
and we finally get: 

w ( U ziOOt ,00 , U ^ A oo t ,00 \ , c s 

Hz = 1 40 00 00 1 00 ^ 2m 2m ) ^ 

which is the only part in Wh that we will keep. Note that there is no contribution from 
1=1, s=l A^ A products. Indeed the Hubbard interaction is invariant under spin rotation 
and couples electrons having zero total S z . As the coefficients of A} A products depend 
solely on / and s they must be zero for s ^ 0. The spectrum for any number of Ti u 
electrons is now easily found: see fig. 2 and table I. 
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V-CONCLUSION 

The ordering of energy levels in the electron-phonon scheme are clearly opposite to 
those of Hund's rule (compare fig.l and fig. 2). The clear signature of what we can call 
"on-ball" pairing is the ground state 1 S of C 6 q\ the two extra electrons are paired by 
the electron-phonon coupling. We note that the U of the Hubbard model appears divided 
by large factors: this is simply due to the fact that the C 60 molecule is large. As a 
consequence, if U ~ 2 eV, Coulomb repulsion may be overwhelmed by phonon exchange. 
With a H g phonon of typical energy 100 meV and coupling 0(1) as suggested by numerous 
calculations 6 ' 7 ' 10 , the quantity A may be tens of meV. 

It seems to us that the cleanest way to probe this intramolecular pairing would be 
to look at solutions of fullerides leading to free anions such as liquid ammonia solutions 
or organic solvents 19-22 . EPR or IR spectroscopy should be able to discriminate between 
the two types of spectra. Measurements by EPR should determine whether or not the two 
extra electrons in CgQ are paired, for example. In near-IR spectroscopy the lowest allowed 
transition for CgQ~ should be at higher energy than that of Cg~ due to the pairing energy 
while in the Coulomb-Hubbard case it is at lower energy. 

Present experiments 19 ' 20 have studied the near-IR spectra of solutions of fulleride 
anions prepared by electrochemical reduction. There are several peaks that do not fit a 
simple Hiickel scheme of levels. They do not have an immediate interpretation in terms of 
vibrational structure 19 ' 20 . With our energy levels in table I, a tentative fit would lead to 
A ps 80 meV assuming U = 0. Such a value leads to intriguing agreement with the major 
peaks seen for CgQ and CgQ while this is no longer the case for CgQ~ and CgQ~. 
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Finally we mention that recent EPR experiements have given some evidence for 
non-Hund behaviour of the fuller ide anions. While one may observe some trends similar 
to the results of the phonon-exchange approximation, it is clear that the model we used 
is very crude. The interplay with conventional Jahn- Teller effect is an important factor 
missing in our study and of a similar order of magnitude. In a bulk conducting solid we 
do not expect the previous scheme to be valid since the levels are broadened into bands. 
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TABLE CAPTIONS: 



Table I: The left column is the electron number. The levels are identified by their quantum 
numbers and the energies are obtained by straightforward perturbation theory. 
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FIGURE CAPTIONS: 



Fig. 1: The levels of fullerene anions taking into account the pho no n- mediated coupling. 



Fig. 2: The energy levels of fullerene anions taking into account a Hubbard interaction. 
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